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The forniulatiou and implementation of ail optimization method called Simultaneous 
Aerodynamic and Structural Design Optimization (SASDO) is shown as an extension of 
the Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) method. It 
is extended by the inclusion of structure element sizing parameters as design variables and 
Finite Element Method (FEM) analysis responses as constraints. The method aims to 
reduce the computational expense incurred in performing shape and sizing optimization 
using state-of-the-art Computational Fluid Dynamics (CFD) flow analysis, FEM struc- 
tural analysis and sensitivity analysis tools. SASDO is applied to a simple, isolated. 3-D 
wing in inviscid flow. Results show that the method finds the same local optimum as 
a conventional optimization method with some reduction in the computational cost and 
without significant modifications to the analysis tools. 


Nomenclature 

b wing semispan 

Cd drag coefficient 

Cl lift coefficient 

C rn pitching moment coefficient 

c r wing root chord 

C( wing tip chord 

F design objective function 

g design constraints 

K stiffness matrix 

L aerodynamic loads 

free-stream Mach number 
n unit normal vector 

p local aerodynamic pressure 

P compliance, the work done by the aerody- 

namic load to deflect the structure 
r/,^ free-stream dynamic pressure 

Q flow-field variables (state variables) at each 

CFD mesh point 
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A< 2 i change in flow solver field variables due to 
better analysis convergence 
AQ2 change in flow solver field variables due to 
design changes 

R aerodynamic state equation residuals at 

each CFD mesh point 

\R/R {) \ norm of the residual ratio, current /initial 
s surface area 

S seniispan wing planform area 

u structural deflections (state variables) 

A t/i change in deflections due to better analysis 
convergence 

A 1/2 change in deflections due to design changes 

W wing weight 

X CFD volume mesh coordinates 

xle location of wing root leading edge 

- chordwise location normalized bv local wing 

c 

section chord 

x t longitudinal location of wing tip trailing 

edge 

z r root section maximum camber 

a free-stream angle-of- attack 

2 design variables 

T structural element size factor 

7 line search parameter 

A operator which indicates a change in a vari- 

able 
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c,c' convergence tolerances for function and gra- 

dients 

A adjoint variable' for Q 

// adjoint variable for u 

r twist angle at wing tip, positive for leading 

edge up 

Subscripts 

d deflected shape 

j jig (undeflected) shape 

■s wing surface mesh 

r wing volume mesh 

Superscripts 

* designates updated value 

gradient with respect to design variables 

Introduction 

Simultaneous Aerodynamic Analysis and Design 
Optimization (SAADO) is a procedure that incorpo- 
rates design improvement within the iteratively solved 
(nonlinear) aerodynamic analysis so as to achieve fully 
converged flow solutions only near an optimal design. 
When SAADO is applied to a flexible wing rather 
than a rigid wing, the linear Finite Element Method 
(FEM) solution is iteratively coupled with the nonlin- 
ear Computational Fluid Dynamics (CFD) solution. 
When design variables that control structural element 
size are included, it is renamed Simultaneous Aerody- 
namic and Structural Design Optimization (SASDO). 
Overall computational efficiency is achieved because 
the many expensive iterative (nonlinear) solutions for 
non-opt imal design parameters are not converged (i.e., 
obtained) at each optimization step. One can ob- 
tain the design iti the equivalent, of a few (rather 
than many) multiples of the computational time for 
a single, fully converged coupled aero-structural anal- 
ysis. SAADO and similar procedures for simultaneous 
analysis and design (SAND) developed by others are 
noted and discussed by Newman et al. 1 These SAND 
procedures appear best suited for applications where 
the discipline analyses involved in the design are non- 
linear and solved iteratively. Generally, convergence 
of these discipline analyses (i.e., state equations) is 
viewed as an equality constraint in an optimization 
problem. From this latter point of view, the SASDO 
method proceeds through infeasible regions of the de- 
sign space which includes not only the design variables 
3, but also the state variables Q and u. A further ad- 
vantage of SASDO is the efficient utilization of existing 
discipline analysis codes (without internal changes), 
augmented with sensitivity or gradient information, 
and yet effectively coupled more tightly than is done in 
conventional gradient- based optimization procedures, 
referred to as nested analysis and design (NAND) 
procedures. 1 A recent overview of aerodynamic shape 


optimization 2 * discusses both NAND and SAND pro- 
cedures in the context of current steady aerodynamic 
opt i m iz at ion research . 

For single-discipline design problems, the distinc- 
tion between NAND and SAND procedures is fairly 
clear and readily seen. With respect to discipline fea- 
sibility (i.e., convergence of the generally nonlinear, 
iteratively solved state equations), these procedures 
can be viewed as accomplishing design by using only 
very well converged discipline solutions (NAND). or 
as converging a sequence of discipline solutions from 
poorly to well as the design progresses (SAND). How- 
ever, the problem formulation and solution algorithms 
may differ considerably. About twenty SAND refer- 
ences are quoted by Newman et al. 1 and Newman et 
ah; 2 these references discuss a variety of formulations, 
algorithms, and results for single-discipline problems 
(mostly CFD applications) in the sense of SAND as de- 
fined above. For mult idisciplinary design optimization 
problems, the distinction between NAND and SAND 
is somewhat blurred because there are feasibility con- 
siderations with respect to all the individual discipline 
state equations, as well as with respect to the mult i- 
disciplinary system compatibility and constraints. A 
number of the papers in Ref. 3 discuss M DO for- 
mulations and algorithms that, are called SAND-like; 
however, not, all of these latter MDO procedures ap- 
pear to agree with the sense of SAND defined above 
and used herein; one that does is Ref. 4. 

Fhe computational feasibility of SAADO for quasi 
1-D nozzle shape design based on the Euler equation 
CFD approximation was demonstrated by Hou et ah 5 
and Mani. h Application of SAADO for turbulent tran- 
sonic airfoil shape design based on a 2-D thin- layer 
Navier- Stokes CFD approximation was demonstrated 
and reported in a later paper by Hou et ah 7 Both of 
these application results are summarized and briefly 
discussed in Ref. 1. The application of SAADO for 
rigid 3-D wing design based on the Euler CFD ap- 
proximation was presented in Ref. 8. These SAADO 
procedures utilized quasi-anaivtical sensitivity deriva- 
tives obtained from hand-differentiated code for the 
initial quasi 1-D application 5 * and from aut omatically 
differentiated code for both the 2-D airfoil applica- 
tion' and the 3-D rigid wing application. 8 Different 
optimization techniques have also been used in these 
SAADO procedures. 

The extension to multidisciplinary analysis with 
shape design variables only was presented in Ref. 9. 
Our initial results from SASDO are given in this pa- 
per. The analysis problem, the objective function, and 
the constraints are the same as those used in Ref. 9. 
That is, changes in design variables are sought to 
produce improvement in the lift-to-drag ratio of a sim- 
ple wing subject to both aerodynamic and structural 
solution-dependent constraints. These constraints are 
the difference between the lift and weight, the pitch- 


2 OF 1 1 

American Institute of Aeronautics and Astronautics Paper 2001-2527 


ing moment coefficient, and the compliance, a function 
representing work done by the aerodynamic load to 
deflect the structure. There are also geometric con- 
straints. 

The flexible wing studied here is formulated as a 
static aeroelastic problem. Similar problems have 
been used as examples in Refs. 10-15 to study var- 
ious solution strategies for multidisciplinary analysis 
and optimization. In particular, Arian 12 analyzed the 
Hessian matrix of the system equations to derive the 
mathematical conditions under which the aeroelastic 
optimization problem can be solved in a “loosely cou- 
pled manner. The multidisciplinary research of Walsh 
et al. 13,14 emphasized engineering aspects of integrat- 
ing high fidelity disciplinary analysis software and dis- 
tributed computing over a network of heterogeneous 
computers. The aeroelastic analysis results of Reuther 
et al. 15 were verified with experimental data. 

Only a limited amount of literature related to aero- 
elastic problems has elaborated on the coupled sen- 
sitivity analysis. Kapania, Eldred and Barthelemy: U) 
Arslan and Carlson; 1 ' and (.Runt a and Sobieszczanski- 
Sobieski 18 derived global sensitivity equations (GSEs); 
some matrix coefficients in these GSEs were evaluated 
by finite differencing. Quinta 1 ' 1 later introduced modal 
coordinates to approximate the elastic displacement 
vector in order to reduce the size of the GSE. Newman, 
Whitfield, and Anderson 20 used the complex variable 
approach to obtain the aeroelastic sensitivity deriva- 
tives, whereas Reuther et al. 15 employed the adjoint 
variable approach to derive the aeroelastic sensitivity 
equations. A mathematical study of the coupled non- 
linear, incompressible aeroelastic analysis and sensitiv- 
ity analysis problems has been given by Gh atlas and 
Li. 21 Recent results on aeroelastic sensitivity analysis 
and optimization can be found in Refs. *2*2-24. Partic- 
ularly, Maute et al. 23 and Hou and Satyanarayana 24 
explicitly formulated the deflection update and the 
load transfer between the separate flow and structures 
solvers as part of the coupled sensitivity equations. 
In the present study, the coupled sensitivity equa- 
tions are constructed by differentiating the aeroelastic 
state equations and solving them by a Generalized 
Gauss- Seidel (GGS) method. 12 The present SASDO 
concept is very similar to t hat of Ghat t as and others, 
Refs. 4,21,25,26, but differs in the implementation 
details as described later. 

Problem Description 

To evaluate the efficacy of the SASDO procedure for 
a problem involving multidisciplinary analysis, it is ap- 
plied herein to a simple, isolated, flexible wing. I he 
wing shape consisted of a trapezoidal planfonn with a 
rounded tip. it was parameterized by fifteen variables; 
five described the planform, and five each described 
the root and tip section shapes. A schematic of the 
wing and its associated shape parameters is shown in 



Fig. 1 Description of semispan wing parameteri- 
zation. 

Fig. 1. The baseline wing section varied linearly from 
an NACA 0012 at t he root to an NACA 0008 at the tip. 
The wing structure consists of a skin, ribs, and spars. 
The ribs and spars consist of shear webs and trusses. 
Six spanwise zones of the structural model are defined 
as depicted in Fig. 1 . The relative sizes of the skin and 
web thicknesses and the truss cross section areas are 
fixed within each zone. Each zone is assigned a param- 
eter T which multiplies all the thicknesses and areas 
of the structural elements in that zone. The specific 
parameters selected as design variables in the sample 
optimization problems are identified in the section en- 
titled Results. The objective function to be minimized 
was the negative of the lift-to-drag ratio, — L/D . Both 
coupled solution-dependent, and geometric constraints 
were imposed. 

The solution-dependent constraints were 

• lower limit on the difference between the total lift 
and the structural weight , Cl * S * q^ — B 

• upper limit on compliance, P = $ pu * nds 

• upper limit on pitching moment, C m , in lieu of a 
trim constraint 

The purely geometric constraints were 

• minimum leading edge radius, in lieu of a manu- 
facturing requirement 

• side constraints (bounds) on the active design 
variables 
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SASDO Procedure 

Formulation 

The flexible SASDO approach formulates the design 
optimization problem as follows: 

min F(Q , Xdv (d< u). ujj) ( 1 ) 

P.Q.u 


subject to 


< 0 : /'= 1,2 i» ( 2 ) 

where the flow field Q and the structural deflection u 
are solutions of the coupled flow equation 

R(Q.XdvU1*u)<ii) = ti (3) 

and the finite element structural equation 

a)) (4) 

The deflected volume mesh, X is determined by 

the deflected surface mesh, X<t 3 , as X ( { v = A’,/, {Xd.s ) • 
I he deflected surface mesh is a result of the jig shape 
augmented by the elastic deflection, u. as Xj? = 
A j(ii) + u. The two disciplines are coupled through 
the deflection, m, and the load, L. 

Recall that Q, /?, and A dv are very large vectors. 
This formulation treats the state variables, Q and u, 
as part of the set of independent design variables, and 
considers the state equations to be constraints. Be- 

cause satisfaction of the equality constraints, Eqs. (3) 
and (4), is required only at the final optimum solution, 
the coupled steady-state aero-structural field equa- 
tions are not converged at every design-optimization 
iteration. The easing of that restriction is expected 
to significantly reduce the excessively large compu- 
tational cost incurred in the conventional approach. 
However, this advantage would likely be offset by the 
very large increase in the number of design variables 
and equality constraint functions, unless some reme- 
dial procedure is adopted. 

Approximations 

I he SASDO method begins with a linearized design 
optimization problem solved for the most favorable 
change in the design variables, Ad, as well as for the 
changes in the state variables, A Q and Au; that is, 


inin 

A i,AQ.Au 


F(Q,AA r (d,w),u,,d) 

+ 5 S- i <?+(Sf+*§fe-)^‘ 

_i_ / PF dX dv \ f i tLL.\ 

' OX \_i, j ' Op ) 

subject to inequality constraints 


Ad 


X/, (d, u), «, d) + 


( 4. a,, 

^ du ^ 0 x, iv ax df ) ** u 


+ 


OQ 


( Mail \ j < n. 

\i*X ttv () A\i s ./ ^ dp J ^ ? 

/ = 1,2 m 


( 5 ) 


( 6 ) 


and equality constraints 


R(Q,Xdv[frv).i3) + %g±Q + 


+ 


( on i 
\o x dv t 




OR OX c 
OX 0 


-Au 


[ 7 ) 


and 


L{Q, X dv {,3 , «))- §kAQ 


K{Xj(i3),t))u 


OQ i 

3 


(*) 


Note that Eqs. (5) through ( 8 ) are linearized approxi- 
mations of Eqs. ( 1 ) through (4), respectively. 

In this formulation, neither the residual of the non- 
linear aerodynamic field equations, B(Q, X , d). nor 
that of the structures equation. K u — L, is required to 
be zero (reach target) until the final optimum design is 
achieved. The linearized problem of Eqs. (5) through 
( 8 ) is difficult to solve directly because of the number 
of design variables and equality constraint equations. 

Direct differentiation method 

One way to overcome this difficulty is by the di- 
rect differentiation method. In this method AQ. Au, 
and Eqs. (7) and ( 8 ) are removed altogether from 
the linearized problem by direct substitution. This 
is achieved by expressing A Q and An as functions of 
AT 

A Q = A Q { + AQ.Ad 

Au = A«i -f AuoA/i ^ ’ 

where vectors AQ\ and Au 1 are corrections in the 
aeroelastic solution due to the improvement of cou- 
pled aeroelastic analysis, while matrices A Qj and Awo 
are corrections due to changes in the design variables. 
These vectors and matrices are solutions of the follow- 
ing coupled sets of equations, obtained from Eqs. (7) 
and ( 8 ): 


on \ r > 


= -/? 


on ox 

0 X dv 0X t i 




A-A«, - + 


10 ) 


where, for the linear FEM, A u — L — 0 at every iter- 
ation, and 




KA,„ = J±A Cl, + jfeflyaj + A.,) Cl) 


( L v' 4 . UK \ u 
{ox, A j + op ) u 


Note that the number of columns of matrices AQ:> 
and Auo is equal to the number of design variables. 
d- Thus the computational cost of Eq. ( 11 ) is directly 
proportional to the number of design variables. 

A new linearized problem with Ad as the only design 
variables can be obtained by substituting Eq. (9) into 
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Eqs. (5) and ( 6 ) for A Q and Aw. 


min 


F(Q, A\ u, i3) + AQi 


dF 

da 


_df^dx A! _\ A 

aArf,, aAj, j -*" 1 


, f wf a ^ i dF dXdv v* 
+ 1 OQ + aA d „ dXT, X :i 


(12) 


+ (*Sfe + £) A »> + *J}^ 


subject to 


<li(Q- A', «, /?) + |^AQi 
+ (*f + 55ftrlfe) i «> 


+ 


{f ^2 + 


^x^ t , aA'd S 


( A j -f A uo 


+ fc i "s + S} A ‘ , S° ; • 


1 , 2 , 


(13) 


The appearance of AQi and Awi in the formulation 
indicates the difference between the SASDO (SAND) 
method and the conventional (NAND) aerodynamic 
optimization method. The AQi and Awi not only 
constitute a change in Q and w, but also play an impor- 
tant role in defining the objective function of Eq. (12) 
and the constraint violation of Eq. (13). We can di- 
rectly solve Eqs. 10 for AQi and Aui, and, in fact, 
in previous SAADO applications that is how AQi was 
determined. However, since AQ^ and Aui, as shown 
in Eq. (10), represent a single Newton's iteration on 
the coupled flow and structures equations, it is possi- 
ble and less computationally expensive to approximate 
their influence on t he solution Q and w by several New- 
ton's iterations of the coupled aeroelastic equations. 
That is, AQi and Aui are not determined explicitly, 
but rather the first three terms of Eqs. 12 and 13 are 
viewed as updated values of F and </, . Note that the 
terms in parentheses in Eqs. (12) and (13) are ap- 
proximated gradients of the objective and constraint 
functions. Once established, this linearized problem 
can be solved using any mathematical programming 
technique for design changes, A j3. Results presented 
in this paper are computed using this direct differen- 
tiation approach. 


Adjoint method 

An alternative way to remove A Q and Aw from the 
linearized problems, Eqs. (5) and ( 6 ). is the adjoint 
variable method. The adjoint variables, A and //, can 
be introduced as the solutions of 


i dR 

1 dQ 


) T A = (f^ 


i>Q 

( dL dXdv \ 

\dx dv dXds J 


F 1 00 


dF \ 
dQ> 


A dX 

( <UL dx dt . \ \ _ ( dF QXdv A __ dF 

^WA'd,, d X d <, ) yWA'rft, WA'd s J du 


dX 

T 


(14) 
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so as to rewrite the expression of the objective function 
in Eq. (5) in terms of A, //, and A 3 as 


F(Q, X 'drill (/.), !/. 3) + A r R + n T (Ku - L) 


+ 
+X T 




dF aAj 


*!■ T 


OR dXdv yt 1 

. dx dl , dx ds ;/ " l_ 


dF \ 

dti ) 

dR \ 

dfi ) 




(15) 

Note that the terms in the brace represent the gra- 
dient of the objective function. The terms A T R and 
fi T (Ku — L) indicate the effect on the design opti- 
mization formulation due to errors in the aeroelastic 
analysis. Furthermore, in the linear sense, the adjoint 
variables and the solution errors can be related by the 
following equat ions: 


dF - r t 
d\- H 


(16) 


and 

~ = (Ku-L) t (17) 

Ofl 

These equations have been mentioned, for example, 
by Pierce and Giles 21 and Venditti and Darmofal 28 
for aerodynamic problems. 

In the typical optimization problem there are many 
design variables. When one can also pose the op- 
timization problem such that there are only a few 
output quantities for objectives and constraints or. in 
the extreme, combine the constraints and objective 
function into a single cost function, the adjoint ap- 
proach to sensitivity analysis has the advantage that 
the adjoint solutions are independent of the number of 
design variables. However, when the disciplines are 
loosely coupled, this approach is impractical since the 
coupled sensitivity analyses would require an adjoint 
for each disciplinary output being transferred, i.e., the 
discretized loads and deflections. In a tightly or im- 
plicitly coupled multidisciplinary analysis the adjoint 
approach may prove practical since this system is anal- 
ogous to a single discipline. 

Line Search 

A one-dimensional search on the step size parame- 
ter 7 is then performed in order to find the updated 
values of A/#, A A', AQ, and Aw. Given the search 
direction A j3 determined by either the direct differen- 
tiation method (Eqs. 12 and 13) or the adjoint method 
equivalent, this line search functions to adjust its mag- 
nitude so as to simultaneously ensure better results for 
both design and analysis (converged solutions). The 
step size parameter 7 plays the role of a relaxation 
factor in the standard Newton's iteration. The search 
procedure employed solves a nonlinear optimization 
problem of the form 


min F(Q* , u (18) 
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subject, to 

<li{Q\X\u\i.r) < 0; ;= 1,2 tn (19) 

R(Q\X\ir) = 0 ( 20 ) 

and 

A* ( A' “ , 3* ) u * = L(Q\ X*) (21) 

where step size 7 is the only design variable. Again 
it is noted for emphasis that the equality constraints. 
Eqs. (20) and (21), are not required to be zero (reach 
target) until the final optimum design; violations of 
these equality constraints must simply be progressively 
reduced until the SASDO procedure converges. 

The updated Q* and u* can be viewed as Q* = Q-\- 
A Q and u* — u + A«. where AQ and Aw satisfy the 
first order approximations to Eqs. (20) and (21). That 
is, AQ and An are the solutions of Eqs. (7) and (8) 
where, in Eq. (9). A3 is replaced by A/?* = 7A/T 
Consequently, Q* = Q + AQ\ + 7 AQ‘ 2 A3 and u* = 
u-hA«i +7 Au*>X3 are readily available once 7 is found. 
The A ( ) 1 terms appearing in the above SASDO for- 
mulation are due to better convergence of the coupled 
analysis, whereas the A()'» terms are due to changes in 
the design variables. In fact, AQ 2 and An *j approach 
the flow field and deflection sensitivities. Q l and u\ as 
the solution becomes better converged. 

Implementation 

The following pseudocode shows algorithmically 
how the method was implemented, 
set initial analysis convergence tolerance. 5 
set initial solution vectors, Q and a, 
set, initial design variables, 3 
do until converged 

1. solve coupled aeroelastic analysis, Eqs. (3) k (4), 
partially converged to z 

2. compute F and g 

3. solve coupled aeroelastic sensitivity analysis. 
Eq. (11). partially converged to s' 

4. compute A3 terms of Eqs. (5) fc (6) 

5. solve optimization problem Eqs. (5) k ((>) for A3 

6. solve Eqs. (18) through (21) for line search pa- 
rameter, 7 

7. update rf, u, and Q 

8. tighten analysis convergence tolerance, 5 = e * 2, 

0 < Z < 1 

end do 

This pseudocode is similar to that used in the 
Biros and Ghattas 25 SAND approach. Specifically, 
both approaches use a Sequential Quadratic Program- 
ming (SQP) method to solve the design equations 
(step 5) and ail approximate factorization method to 
solve the system equations (step 1). Step 3 above uses 
an incremental iterative method with approximate 
factorization to solve for derivatives in direct mode 
rather than as a solution of the adjoint equation of 



Fig. 2 Diagram of SASDO procedure. 

Biros and Ghattas. Jo In addition, the line search step 
(step 6) and the convergence tightening step (step 8) 
were not included in the Biros and Ghattas method. 
A schematic of the present SASDO procedure is shown 
in Fig. 2. The dashed box, labeled "Partially Con- 
verged System Analyis,” depicts the coupled analysis 
iteration loop. Steps 1 and 2 of the pseudocode; that 
labeled "Partially Converged Sensitivity Analysis' 
depicts the coupled derivative iteration loop. Step 3; 
that labeled “Partially C onverged Design” depicts the 
design steps, Steps 0-8 of the pseudocode. Specific 
computational tools and methods used to perform 
the tasks depicted by the solid boxes in Fig. 2 are 
identified in the next section. 

Computational Tools and Models 

Major computations in this SASDO procedure are 
performed using a collection of existing codes. These 
codes are executed by a separate driver code and 
scripts that implement the SASDO procedure as just 
discussed. Each code runs independently, some simul- 
taneously, on different processors, and the required 
I/O transfers between them, also directed by the 
driver, are accomplished via data files. 

The aerodynamic flow analysis code used for this 
study is a version of the OFL3D code. 3,1 Only Eu- 
ler analyses are performed for this work, although the 
code is capable of solving the Navier-Stokes equations 
with any of several t urbulence models. The gradient 
version of this code, which was used for aerodynamic 
sensitivity analysis, was generated by an unconven- 
tional application 31 of the automatic differentiation 
code ADIFOR 32 33 to produce a relatively efficient, 
direct mode, gradient analysis code. CFL3D.ADII. 34 
It should be pointed out that the ADIFOR. process 
produces a discretized derivative code consistent with 
the discretized function analysis code. The addition of 
a stopping criterion based on the norm of the residual 
of the field equations was the only modification of the 
CFL3D.ADII code made to accommodate the SASDO 
procedure. 
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CFD mesh 
C-O topology 
73x25x25 volume 
49x25 on the wing 

FEM mesh 
3251 elements: 

2240 CST 

1011 truss 
583 points 

Fig. 3 Computational meshes for wing analysis. 

The surface geometry was generated based on the 
parameters described in a previous section by a code 
utilizing the Rapid Aircraft Parameterization Input 
Design (RAPID) technique developed by Smith, et. 
al. 35 This code was preprocessed with ADIFOR to gen- 
erate a code capable of producing sensitivity deriva- 
tives as well. 

The CFD volume mesh needed by the flow analysis 
code was generated using a version of the CSCMDO 30 
grid generation code. The associated grid sensitivity 
derivatives needed by the flow sensitivity analysis were 
generated with an automatically differentiated ver- 
sion of CSC MIX). 3 ' In addition to the parameterized 
surface mesh and accompanying gradients. CSCMDO 
requires a baseline volume mesh of similar shape and 
identical topology. The 45,000 grid point baseline vol- 
ume mesh of C-0 topology used in the present flexible 
wing examples was obtained with the Gridgen™ code. 
The wing surface portion of the mesh is shown in 
Fig. 3. This mesh is admittedly particularly coarse 
by current CFD analysis standards. 

The structural analysis code 38 used to compute 
the deflection of the elastic wing was a generic fi- 
nite element code. The flexible structure for the wing 
shown in Fig. 3 was discretized by 583 nodes; there 
were 2,141 constant -strain triangle (CST) elements 
and 1,110 truss elements. Zone boundaries for the de- 
sign variables controlling element size are also shown in 
Fig. 3. Because the elastic deformation was assumed 
to be small, linear elasticity was deemed to be ap- 
propriate. The structural sensitivity equations were 
derived based upon the direct differentiation method. 
We note that the sensitivity of the aerodynamic forces 
appears as a term on the right-hand side (RHS) of the 
deflection sensitivity equations. The derivative of the 
stiffness matrix in these sensitivity equations was also 
generated 3 ^ by using the ADIFOR 33 ' 33 technique. We 
note that the coefficient matrix of the structural sensi- 
tivity equations was identical to that of the structural 
equations. Consequently, these structural sensitivity 
equations were solved efficiently by backward substi- 
tution with different RHSs for each sensitivity. 

At the wing surface, i.e., the interface where aero- 
dynamic load and structural deflection information is 
transferred, it was assumed that surface nodes of the 



FEM structural model were a subset of the CFD aero- 
dynamic surface mesh points (see Fig. 3) for the 
present SASDO application. This lack of generality 
allowed for simplifications in the data transfers and, al- 
though an important issue, it was not deemed crucial 
for these initial SASDO demonstrations. Future ap- 
plications to more complex configurations should allow 
for transfer of conserved information between arbitrary 
meshes as required by the individual disciplines. A 
recent review of such data transfer techniques and a 
specific proposed one are given in Ref. 40. 

Conventional (NAND) and SASDO (SAND) proce- 
dures were implemented using the SQP method of t he 
DOT 11 optimization software. All computations were 
executed on an SCR Origin 2000™ workstation with 
250MHz R 10000™ processors. The CFD sensitiv- 
ity calculations were partitioned and run on several 
processors to reduce required memory and elapsed op- 
timization time. This partitioning, however, results in 
additional accumulated computation time due to the 
nat ure of ADIFOR-generated sensit ivity analysis code. 

Results 

The optimization results shown in this work are 
for design problems involving only four or eight de- 
sign variables out of the 21 available wing parameters. 
The results shown by these authors in an earlier work 51 
used design variables that directly affected either the 
aerodynamic analysis alone or both the aerodynamic 
analysis and the structural analysis. In this work, ad- 
ditional design variables are chosen that directly affect 
only the structural analysis. The flow conditions for 
the wing optimizations were — 0.8 and o = 1°. 



Fig. 4 Comparison of planforni shapes and sur- 
face pressure contours for 4-design- variable cases. 

Moo = 0.8, a = 1°. 

Four-Design-Variable Problems 

Table 1 and Fig. 4 show results of several optimiza- 
tion problems involving four design variables: the tip 
chord c f , the tip setback jv, and the structural element 
size factor for the two most inboard zones, \\ and 
I 2 . Two of the cases, designated Convc and SASDOe 
in Table 1, represent direct comparisons of SASDO 
and the conventional method for consistent accuracy 
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Table 1 Summary of Four-Design- Variable Cases. 



baseline 

SASDOa 

Optimizations 
Convn SASDOc 

Convc 

C T 

L 

1.07778 

1.07281 

1.07827 

1.07697 

Xt 

1. 

1.93259 

1.91880 

1,93052 

1.92614 

I’l 

1. 

1.06631 

1.08008 

1.09637 

1 .09982 

I'o 

i. 

0.69510 

0.69187 

0.68664 

0.68691 

F 

-7.149 

-10.187 

-10.181 

-10.186 

-10.183 

gi 

-.0302 

-.000760 

-.00177 

-.0000629 

-.000915 

g2 

-.8882 

-0.880 

-0.876 

-0.880 

-0.879 

g3 

-.2647 

-.000639 

-.000159 

-.000149 

+.000108 

c 

9.6e-7 

9.9e-5 

8.9e-7 

6.6e-6 

4.0e-6 

Y state 

1 

1.1 

4.8 

2.0 

3.9 

Es rad 

— 

27.1 

63.1 

26.9 

40.4 


of funct ion and gradient analyses. The other t wo cases 
show effect s of changes in the accuracy of the function 
and gradient analyses. The resulting designs are es- 
sentially identical for all four cases. Figure 4 shows 
a comparison of the wing planform and the surface 
pressure coefficient results for the baseline design and 
the design designated SASDOf in Table 1. The shock 
wave has been weakened somewhat in the optimized 
cases from that on the original wing, as one would 
expect. As one can see qualitatively in Fig. 4 and 
numerically from the values of the objective function 
/’, the constraints </;, and the final design variables 
in fable 1, the final designs are very similar for the 
four problems. The relative computational costs of 
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> 
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JJ 

id 

rH 

© 



conven- 

tional 


SASDO 



d(mesh) 

d(FEM) 

d(CFD) 

mesh 

FEM 

CFD 


Fig. 5 Comparison of computation cost of four- 
design-variable optimization problem using the 
conventional and SASDO methods. 


the optimizations are shown in Table 1 and Fig. 5. 
The accumulated function and gradient analysis times 
(denoted as Estate and S grad ) are shown separately. 
They have been normalized by the cost of the base- 
line coupled function analysis. The components of the 
analyses were lumped together in Table 1 but they 


have been separated in Fig. 5. The total time for 
performing this optimization problem was reduced by 
36 percent using the SASDO method. The analysis 
alone was reduced by 55 percent, but the gradient eval- 
uation was the dominant cost. 



Fig. 6 Comparison of planform shapes and sur- 
face pressure contours for 8-design-variable cases, 
A/.x = 0.8, o = l°. 


Eight-Design- Variable Problems 

Table 2 and Fig. 6 show results of three opt imization 
problems involving eight design variables: the same set 
used in the four-design- variable cases with the inclu- 
sion of the span 6, the root section max camber z r , 
and the structural element size factor for two more 
zones, T3 and T 4. Two of the cases, designated Coiivb 
and SASDOb in Table 2, represent direct comparisons 
of SASDO and the conventional method for consis- 
tent accuracy of the function and gradient analyses. 
Figure 6 shows a comparison of wing planform and 
surface pressure coefficient results for the baseline de- 
sign and the design designated SASDOb in Table 2. 
The relative computational costs of the optimizations 
are shown in Table 2 and Fig. 7. The total time for 
performing this optimization problem was reduced by 
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Table 2 Summary of Eight-Design- Variable Cases. 



baseline 

Optimizations 

Conv A SASDOb Convn 

Ct 

1 . 

0.76574 

0.75929 

0.75794 

Xt 

1 . 

2.(30000 

2.60000 

2.60000 

b 

1 . 

1.05992 

1 .0754:$ 

1.07269 

z r 

1 . 

1.11511 

1.11173 

1.11175 

r i 

1 . 

2.42444 

2.55263 

2.44498 

r 2 

1 . 

2.09537 

2.09924 

2.09394 

r 3 

1. 

1.03(335 

1.01883 

1.05778 

r 4 

1 . 

0.(35137 

0.81046 

0.76612 

F 

-7.149 

-21.255 

-21.235 

-21.244 

gi 

-.0302 

-.7536 

-.7255 

-.7316 

g‘2 

-.8882 

-.00336 

-.00671 

-.000241 

g3 

-.2647 

-.000420 

-.000962 

-.000134 

cr 

9.6e-7 

6.8e-7 

1.8e-5 

U6e-5 

Y state 

1 

13.5 

4.4 

11.3 

£ s rad 

- 

146.7 

1 15.9 

148.6 



conven- SASDO 
tional 


Fig. 7 Comparison of computation cost of eight- 
design-variable optimization problem using the 
conventional and SASDO methods. 

26 percent using the SASDO method. The analysis 
alone was reduced by 60 percent, but, as with the four- 
design- variable problem, the gradient evaluation was 
the dominant cost. 

Further Discussion 

The relative costs, based on CPU timing ratios, 
for the SASDO (SAND) procedures applied to these 
present small 3-D aerodynamic /structural design opti- 
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mization problems are about seven-tent hs of the costs 
of the corresponding conventional (NAND)procedures. 
This range is very similar to that reported for 2-D 
nonlinear aerodynamic shape design optimization in 
Refs. 1 and 4, even though many of the computational 
details differ. The results given in Ref. 1 were for a 
turbulent transonic flow with shock waves computed 
using a Navier-Stokes code; a direct differentiation ap- 
proach (using ADI FOR ) was used for the sensitivity 
analysis. The results reported in Ref. 4 were for a com- 
pressible flow without shock waves computed using a 
nonlinear potential flow code; an adjoint approach was 
used for the sensitivity analysis. Since these two opti- 
mization problems were also not the same, no timing 
comparison between these adjoint and direct, differ- 
entiation solution approaches would be meaningful. 
As indicated earlier, an expected speed-up for using 
an adjoint approach instead of the direct, differentia- 
tion approach was estimated in Ref. 1. (Riattas and 
Bark 20 recently reported 2-D and 3-D results for op- 
timal control of steady incompressible Navier-Stokes 
flow that demonstrate an order-of-magnitude reduc- 
tion of CPU time for a SAND approach versus a 
NANI) approach. These results were obtained using 
reduced Hessian SQP methods that avoid converging 
the flow equations at each optimization iteration. The 
relationship of these met hods wit h respect, to other o|>- 
timization techniques is also discussed in Ref. 26. The 
“Control Theory” approach of Jameson 2 *’ and several 
other SAND-like methods for simultaneous analysis 
and design, which were summarized and discussed by 
Ta'asan. 4 " have been applied to aerodynamic shape 
design problems at several fidelities of CFD approx- 
imation. These techniques have obtained an aerody- 
namic design in the equivalent of several analysis CPU 
t imes for some sample problems. 

Concluding Remarks 

This study has introduced an implementation of the 
SASDO technique for a simple, isolated wing. Initial 
results indicate that SASDO 

1. is feasible under dual simultaneity (i.e. simultane- 
ity not only with respect to analysis and design 
optimization, but also simultaneity with respect 
to flexible wing aero-structural interaction) 

2. finds the same local minimum as a conventional 
technique 

3. is computationally more efficient than a conven- 
tional gradient, -based optimization technique 

4. requires few modifications to the analysis and sen- 
sitivity analysis codes involved 

5. is effective at reducing the function analysis cost, 
but the gradient analysis time is the dominant 
cost 
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